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In heavy ion collisions, event-by-event fluctuations in participating nucleon positions can lead 
to triangular flow. With fluctuating initial conditions, flow coefficients will also fluctuate. In a 
hydrodynamic model, we study the fluctuations in elliptic and triangular flow, due to fluctuating 
initial conditions. Both elliptic and triangular flow fluctuates strongly, triangular flow more strongly 
than the elliptic flow. Strong fluctuations greatly reduce the sensitivity of elliptic and triangular 
flow to viscosity. 

PACS numbers: 47.75+f, 25.75.-q, 25.75.Ld 



I. INTRODUCTION 

It is expected that collisions between two nuclei at 
ultra-relativistic energies will lead to a phase transition 
from hadrons to the fundamental constituents, quarks 
and gluons, usually referred to as Quark-Gluon-Plasma 
(QGP). One of the experimental observables of QGP is 
the azimuthal distribution of produced particles. In non- 
zero impact parameter collision between two identical 
nuclei, the collision zone is asymmetric. Multiple col- 
lisions transform the initial asymmetry into momentum 
anisotropy. Momentum anisotropy is best studied by de- 
composing it in a Fourier series, 
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<j) is the azimuthal angle of the detected particle and tp n 
is the plane of the symmetry of initial collision zone. 
For smooth initial matter distribution, plane of sym- 
metry of the collision zone coincides with the reaction 
plane (the plane containing the impact parameter and 
the beam axis), i/) n = ^ap,Vn. The odd Fourier co- 
efficients are zero by symmetry. However, fluctuations 
in the positions of the participating nucleons can lead 
to non-smooth density distribution, which will fluctuate 
on event-by-event basis. The participating nucleons then 
determine the symmetry plane (tppp), which fluctuate 
around the reaction plane As a result odd harmon- 
ics, which were exactly zero for smoothed initial distri- 
bution, can be developed. It has been conjectured that 
third harmonic V3, which is response of the initial tri- 
angularity of the medium, is responsible for the observed 
structures in two particle correlation in Au+Au collisions 
0,0, 0|, [|, 0,0. The ridge structure in pp collisions 
also has a natural explanation if odd harmonic flow de- 
velops. Recently, ALICE collaboration has observed odd 
harmonic flows in Pb+Pb collisions |8|. In most central 
collisions, the elliptic flow («a) and triangular flow (v^) 
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are of similar magnitude. In peripheral collisions how- 
ever, elliptic flow dominates. 

The second harmonic or the elliptic flow (i^) has been 
studied extensively in y/s NN =200 GeV Au+Au collisions 
at RHIC 0, GlJ. Recently, ALICE collaboration mea- 
sured elliptic flow in ^^^=2.76 TeV Pb+Pb collisions 
at LHC [1, [Tll| - Large elliptic flow has provided com- 
pelling evidence that at RHIC and LHC, nearly perfect 
fluid is produced. Deviation from the ideal fluid be- 
havior is controlled by shear viscosity to entropy ratio 
(t)/s). Effect of shear viscosity is to dampen the flow 
coefficients. Elliptic flow has sensitive dependence on 
77/s. In smooth hydrodynamics, sensitivity of elliptic flow 
has been utilised to obtain phenomenological estimates 
of 77/s [I3 - IT9j |. Triangular flow is supposed to be more 
sensitive to viscosity than the elliptic flow 0, 0] and one 
expects that triangular flow measurements will constrain 
77/s more accurately. 

Evcnt-by-event fluctuations in initial conditions gen- 
erate the triangular flow. It is then natural that the 
triangular flow itself will also fluctuate, event-by-event. 
Unless the fluctuations are within some reasonable limit, 
sensitivity of the flow to rj/s will reduce greatly. In the 
present paper, in a hydrodynamic model, we have stud- 
ied the fluctuations in elliptic and triangular flow due to 
fluctuating initial conditions. Its sensitivity to the 77/s is 
also studied. It appear that with fluctuating initial con- 
ditions, the sensitivity of elliptic and triangular flow to 
77/s is greatly reduced. Triangular flow fluctuates more 
strongly than the elliptic flow and become even less sen- 
sitive to 77/s than the elliptic flow. Large fluctuations 
belie the possibility of constraining viscosity to entropy 
ratio from triangular flow measurements. 



II. HYDRODYNAMIC EQUATIONS, 
EQUATION OF STATE AND INITIAL 
CONDITIONS 

We assume that in -^^^=2.76 TeV, Pb+Pb collisions 
at LHC, a baryon free fluid is formed. Only dissipative 
effect we consider is the shear viscosity. Heat conduction 
and bulk viscosity is neglected. The space-time evolution 
of the fluid is obtained by solving, 
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FIG. 1: Distribution of spatial eccentricity (62) and triangu- 
larity (£3) in b=8.9 fm Pb+Pb collisions in a model with hot 
spot formation. Hot spot width a=l fm. 
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FIG. 2: comparison of elliptic and triangular flow in hydro- 
dynamical simulations with -/V eue „t=2500, 500 and 50. Data 
points are marginally shifted to distinguish between them. 



d^" = 0, (2) 
Dtt^ = - — (7^ - 2»/V <, V > ) 

- [tAr"* + u v ti^\Du x . (3) 

Eq[2] is the conservation equation for the energy- 
momentum tensor, T^ v = (e + p)u^u v — pg^ v + n^ v , 
e, p and u being the energy density, pressure and fluid 
velocity respectively. ■k^ v is the shear stress tensor. EqJ3] 
is the relaxation equation for the shear stress tensor tt^ v . 
In EqJ2l D — u /1 9„ is the convective time derivative, 
V<»u v> = i(Vu" + Vu^) - \{d.u){g> lv - uV) is a 
symmetric traceless tensor. r\ is the shear viscosity and t„ 
is the relaxation time. It may be mentioned that in a con- 
formally symmetric fluid relaxation equation can contain 
additional terms [§(|. Assuming boost- invariance, Eqs|3] 
and [3] are solved in (r = y/t 2 — z 2 , x, y, i] s = glnfzf) 
coordinates, with the code " 'AZHYDRO-KOLKATA'", 
developed at the Cyclotron Centre, Kolkata. Details of 
the code can be found in [21] . 

Eqs l2l3l are closed with an equation of state (EoS) p = 
p(s). Currently, there is consensus that the confinement- 
deconfincmcnt transition is a cross over and the cross 
over or the pseudo critical temperature for the transi- 
tion is T c fa 170 MeV [HHIIj]. In the present study, we 
use an equation of state where the Wuppertal-Budapest 
[22I Izij lattice simulations for the deconfined phase is 
smoothly joined at T = T c = 174 MeV, with hadronic 
resonance gas EoS comprising all the resonances below 
mass m res =2.5 GeV. Details of the EoS can be found in 

Solution of partial differential equations fEqs l2!3p re- 
quires initial conditions, e.g. transverse profile of the 
energy density (e(x,y)), fluid velocity (v x (x,y),v y (x,y)) 
and shear stress tensor (it 1 - 1 " (x,y)) at the initial time 
Tj. One also need to specify the viscosity (17) and the 
relaxation time (t„). A freeze-out prescription is also 
needed to convert the information about fluid energy 
density and velocity to particle spectra. We assume 



that the fluid is thermalised at Ti=0.6 fm and the ini- 
tial fluid velocity is zero, v x (x,y) = v y (x,y) — 0. We 
initialise the shear stress tensor to boost-invariant val- 
ues, ir xx — ir vv — 2rj/3ri, tt xv —0 and for the relaxation 
time, we use the Boltzmann estimate tv = 3r)/2p. We 
also assume that the viscosity to entropy density (77/s) 
remains a constant throughout the evolution and simu- 
late Pb+Pb collisions for a range of 77/s. The freeze-out 
is fixed at T F =130 MeV. 

For the fluctuating initial energy density distribution, 
we use a model of hot spots in the initial states. Re- 
cently, similar models are used to study elliptic flow in 
pp collisions at LHC [§6-28]. In [2(|, a similar model 
was used to study anisotropy in heavy ion collisions due 
to fluctuating initial conditions. We assume that in an 
impact parameter b collision, each participating nucleon 
pair randomly deposit some energy in the reaction vol- 
ume, and produces a hot spot. The hot spots are assumed 
to be Gaussian distributed. The initial energy density is 
then super position of N = N part icipant hot spots. 

^participant ^ ^ 

e(x,y)=e e (4) 

i=l 

The participant number N part i C ipant is calculated in a 
Glauber model. We also restrict the centre of hotspots 
(fi) within the transverse area defined by the Glauber 
model of participant distribution. The central density 
£0 and the width a are parameters of the model. We 
fix er=l fm. The central density £0 is fixed to reproduce 
approximately the experimental charged particles in a 
peripheral (30-40%) Pb+Pb collisions. 

Since hot spot positions are random, density distribu- 
tion fluctuates from event to event. Asymmetry in the 
initial energy density distribution can be characterised in 
terms of e n and ip n [j| , 

ene in^ n = J J e(x,y)r 2 e m *dxdy 
J J s(x, y)r 2 dxdy 
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FIG. 3: (color online) Integrated elliptic and triangular flow 
for viscosity to entropy ratio r//s=0, 0.04, 0.08, 0.12 and 0.16. 
Event size is N event = 500. 



Here, e„'s are parameters characterizing the asymme- 
try in the distribution, e.g. dipole asymmetry: e\, partic- 
ipant eccentricity: £2, participant triangularity: £3 etc. 
tp n is the polar angle of the participant plane for n-th 
order flow. Teaney and Yan Q, from theoretical consid- 
eration, argued that for triangular flow, the r 2 term in 



EqJS] should be replaced by 



£3e 



/ / e(x,y)r 3 e l3( ^dxdy 
J J e(x, y)r 3 dxdy 



(6) 



In the present simulation, we have used EqJBJto com- 
pute triangular flow. Probability distribution of spa- 
tial eccentricity (£2) and triangularity (£3), in b=8.9 fm 
Pb+Pb collision is shown in Fig[TJ Eccentricity distribu- 
tion peaks around £2 ~ 0.2. The triangularity parameter 
£3 peaks around £3 « 0.04. The fluctuations of £3 however 
is much larger than that in £ 2 . Since elliptic and trian- 
gular flows are response respectively to the initial state 
eccentricity and triangularity, larger fluctuations are ex- 
pected in triangular flow than in elliptic flow. 



III. VISCOUS EFFECTS ON ELLIPTIC AND 
TRIANGULAR FLOW 

For fluid viscosity to entropy ratio rj/s=0, 0.04, 0.08, 
0.12 and 0.16, we have simulated b=8.9 fm Pb+Pb col- 
lisions. b=8.9 fm collisions approximately corresponds 
to 30-40% collision. In viscous evolution, entropy is 
generated. To account for the entropy generation, the 
Gaussian density e$ was reduced with increasing viscos- 
ity, such that in ideal and viscous fluid, on the average, 
7r~ multiplicity remains the same. In each event, Israel- 
Stewart's hydrodynamic equations are solved and from 
the freeze-out surface, invariant distribution ( , d S — ) for 

7r~ was obtained. In analogy to EqjSl invariant distribu- 
tion can be characterised by 'harmonic flow coefficients' 

M. 
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FIG. 4: (color online) Differential elliptic flow (1)2) for viscos- 
ity to entropy ratio r)/s=0, 0.04, 0.08, 0.12 and 0.16. Event 
size is iVeuent =60. Data points are shifted marginally to dis- 
tinguish between them. 
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a intj) dN 
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dN 



(7) 



/ P Td PT d(j)e in * ^ 



dypT dpT dtp 



<IN 
dy 



(8) 



In a boost-invariant version of hydrodynamics, flow 
coefficients are rapidity independent and in the follow- 
ing, we have dropped the rapidity dependence. Present 
simulations are suitable only for central rapidity, y ~0, 
where boost-invariance is most justified. In the present 
study, we have used N event =500 events. Is sample size 
N event =b§§ sufficiently large to comment conclusively on 
fluctuations in w 2 and V3? In FigJ^l simulation results for 
the elliptic (t^) and triangular (v^) flow are shown. The 
filled (open) circles, squares and triangle corresponds to 
elliptic (triangular) flow for event size N event =2500, 500 
and 50 respectively. The fluid is assumed to be ideal. 
The symbols represent event average and the error bars 
the variance. If fluctuations are of statistical origin, one 
expect the variance in vi and v 3 increase by a factor of ~ 
2 between event size 2500 and 500 and by a factor of ~7 
between event size 2500 and 50. No such increase is seen 
in simulations. As shown in Fig(5]the event averaged el- 
liptic and triangular flow are approximately independent 
of the event size. More importantly, the variances are 
also approximately independent of event size. The result 
confirms that fluctuations are systematic rather than sta- 
tistical. 

We may mention here that the present model of fluctu- 
ating initial conditions do not induce large fluctuations in 
the 7T~ transverse momentum spectra, px spectra remain 
largely unaffected. The model is also consistent with the 
experimental observation that in central collisions, ellip- 
tic and triangular flows are of similar magnitude, but in 
peripheral collisions, elliptic flow dominates. 

In Fig [31 viscosity dependence of integrated elliptic and 
triangular flow is shown. With initial conditions both 
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FIG. 5: (color online) same as in Fig[4]but for triangular flow. 



V2 and V3 fluctuates. However, fluctuations in V2 is not 
large, ~20%. Within the fluctuations, integrated vi does 
not show any viscosity dependence. Triangular flow fluc- 
tuates more strongly than the elliptic flow. Uncertainty 
in t>3 is ^50%. Larger fluctuations in is not unex- 
pected. As shown in FigQ] triangularity fluctuates more 
strongly than the eccentricity. It is then expected that 
the will fluctuate more strongly than v-i- Integrated 
triangular flow also does not show any viscosity depen- 
dence. 

In FigH]and[51 differential elliptic and triangular flow is 
shown as a function of viscosity. Average value of V2(j>t) 
and v 3 (pt) decreases with increasing viscosity, however, 
both elliptic and triangular flow fluctuates strongly. In 
the pt range 1-3 GeV, fluctuations in V2{pt) is ~15-20%. 
Fluctuations in v%{pt) is even more ~ 70-80%. The fluc- 
tuations in anisotropic flow greatly reduce their efficacy 
as a diagnostic tool. For example, within the fluctua- 
tions, differential elliptic flow does not distinguish be- 
tween ideal fluid and fluid with viscosity to entropy ratio 
0.08. Triangular flow is even more insensitive. Fluid vis- 
cosity varying between 0-0.16 is not distinguished. 

Triangular flow depends on the definition of triangu- 
larity. In the above simulations, we have used the defini- 
tion due to Teaney and Yan 0]. As noted earlier, Alver 
and Rolland [5[ used an alternate definition for trian- 
gularity (see EqO. It is interesting to compare v 3 in 
the two definitions. In Fig[5J for fluid viscosity ?y/s=0.08 
and 0.16, triangular flow from the two definitions is com- 
pared. Average triangular flow and its fluctuations are 
approximately similar in both the definitions. The re- 
sult is understood. Triangular flow depend explicitly on 
the participant plane angle ^3. 1P3 is not changed much 
between the two definition. For example, in the present 
simulation, out of 500 events, in ~400 events ipz in two 
definitions differ by less than 10%. Marginal change in 
ips lead to similar triangular flow. 

Can the fluctuations in ^3 be reduced? As demon- 



strated earlier, larger event size do not reduce the fluctu- 
ations. One of the parameter of the model is the Gaus- 
sian width a of the hot spots. Initial density distribution 
become more asymmetric if a is reduced. We have used 
cr=l fm. Explicit simulations indicate that if a is reduced 
to 0.5 fm, in the px range 1-3 GeV, average of triangular 
flow increase by a factor of ^2, the fluctuations however 
remains largely unaltered. In a more realistic simulation, 
impact parameter would also fluctuate and fluctuations 
in V3 will be even more. Indeed, if the initial state fluc- 
tuations are the 'only' source of the triangular flow, large 
fluctuations in the triangular flow is inevitable. Fluctu- 
ations will reduce if, apart from initial triangularity, jet 
quenching, Cerenkov radiation etc. also contribute to the 
triangular flow. Indeed, hydrodynamical simulations in- 
dicate that triangularity in density distribution can be 
developed due to jet quenching (see Fig. 2 of ref.[3fj), 
whose response would be the triangular flow. 
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FIG. 6: (color online) Comparison of triangular flow with 
definition of participant plane due to Teaney and Yan [?J and 
Alver and Rolland [jj. Simulation results for fluid viscosity 
over entropy ratio, ?7/s=0.08 and 0.16, are shown. 



IV. CONCLUSIONS 

In a simple model of fluctuating initial states, we have 
studied fluctuations in elliptic and triangular flow in ideal 
and viscous fluid evolution. It is shown that fluctuations 
in triangular flow can be very large. Large fluctuation re- 
duces the sensitivity of flow coefficients to viscosity. For 
example, fluctuations of elliptic flow make it insensitive 
to variation of viscosity to entropy ratio in the range 
7//s=0-0.08. Fluctuations in triangular flow are even 
larger and viscosity to entropy ratio varying between 0- 
0.16 is not distinguished. We conclude that if the initial 
state fluctuations are the only source of triangular flow, 
triangular flow will be greatly insensitive to viscosity. 
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